Finite size analysis of a two-dimensional Ising model within a nonextensive approach 
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In this work we present a thorough analysis ol the phase transitions that occur in a ferromagnetic 2D 
Ising model, with only nearest-neighbors interactions, in the framework of the Tsallis nonextensive 
statistics. We performed Monte Carlo simulations on square lattices with linear sizes L ranging 
from 32 up to 512. The statistical weight of the Metropolis algorithm was changed according to 
the nonextensive statistics. Discontinuities in the m(T) curve are observed for q < 0.5. However, 
we have verified only one peak on the energy histograms at the critical temperatures, indicating 
the occurrence of continuous phase transitions. For the 0.5 < q < 1.0 regime, we have found 
continuous phase transitions between the ordered and the disordered phases, and determined the 
critical exponents via finite-size scaling. We verified that the critical exponents a, /3 and 7 depend 
on the entropic index q in the range 0.5 < q < 1.0 in the form a{q) = (10 5^ — 33g -I- 23)/20, 
P{q) = {2q — l)/8 and 7(5) = (g^ — q + 7)74. On the other hand, the critical exponent v does not 
depend on q. This suggests a violation of the scaling relations 2/3-1-7 = and a-\-2/3-\-j = 2 
and a nonuniversality of the critical exponents along the ferro-paramagnetic frontier. 

PACS numbers: 05.10.Ln, 05.50.+q, 05.70.Fh, 05.90.-|-m, 



I. INTRODUCTION 



Inspired in the geometrical theory of multifractals. Tsalhs has suggested a generalization of the Boltzmann-Gibbs 
entropy (Sbg), which is known as the nonadditive entropy P, The entropy form is postulated to be 

s, = , (1) 

where '^^Pi — 1 and fc is a constant. The idea behind this generalization is that Sq is the measure of the information 
of biased systems. Thus, being pi the probability of finding a given system on the state i, the factor q introduces a bias 
into the probabihty set, i.e., if < < 1 then pj > pi for q < 1 and pf < pi for q > 1. In other words, q < I privileges 
the less probable events in opposition to the more probable ones, and vice-versa. This entropy is invariant under 
permutations, becomes zero for the maximum knowledge about the system, and for g = 1, i.e. for unbiased systems, 
it would recover the BG entropy. The bias factor q is called the entropic index and g £ 5R. Recently, it has been 
proposed that q is connected to the dynamics of the system 0, 0, HI, S 0, [1] ■ Besides representing a generalization, 
the nonextensive entropy 5*^, as much as Sbg, is positive, concave and Lesche-stable {\/q > 0). It has also been shown 
that for systems with certain types of correlations that induces scale invariance in the phase space, the entropy Sq 
becomes additive [13, [O, [13 • The optimization of the entropy in Eq. ([J) leads to the equilibrium distribution and a 
generalization of the Boltzmann-Gibbs statistics, that is called nonextensive statistics. This generalization has been 
successfully applied in many areas of physics, biology and computation in the past few years [l3l.[l^. [Tsl. [iGj. 

On the other hand, magnetic models are one of the most studied systems in condensed matter, and the Ising model 
is a prototype that has been extensively investigated for the last 30 years. More recently, some works have investigated 
the magnetic properties of some manganese oxides, called manganites, and connections with the nonextensive statistics 
have been proposed [1, [O, [H, [ill ■ ^ recent work, Rcis et al. [ll] studied the phase transitions that occur in a 
classical spin system within the mean- field approximation, in the framework of Tsallis nonextensive statistics, and some 
interesting properties were found: the system presents first-order phase transitions for q < 0.5, but only continuous 
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transitions were found for 0.5 < q < 1.0. The results present qualitative agreement with experimental data in the 
Lao.eoYo.orCao.aaMnOa manganite [20 |. 

From this, two natural questions appear: What properties of infinitc-rangc-intcraction models can appear in short- 
range-interaction ones? Can the mean-field predictions be verified on low dimensional systems? Although a 2D 
Ising model, defined in the limit of nearest-neighbors interactions, is usually treated according the BG statistics, a 
nonextensive approach can be seen as a toy model for the elucidation of the question: will the mean-field behavior 
be also present in a short-range-interaction model? Furthermore, using this kind of study one can also verify the 
accuracy of the nonextensive statistics applied to magnetic systems. 

In a attempt to clarify some of these questions, we report in this paper results on the study of the phase transitions 
that occur on a 2D Ising model within a nonextensive approach. These results were obtained through Monte Carlo 
(MC) simulations upon replacing the statistical weight of the Metropolis algorithm by the nonextensive one. We 
performed a finite size scaling in order to estimate the critical exponents for different values of the entropic index 
q S [0, 1]. As discussed in Ref . [2]| . the critical temperatures Tc depend on g (Vq G [0, 1]), but we have found in the 
present work that the critical exponents a, (3 and 7 depend on the entropic index q, in the range 0.5 < q < 1.0. 



II. NONEXTENSIVE STATISTICS AND MONTE CARLO SIMULATION 

In nonextensive statistics theory (see, e.g., Refs.p. [22l. [23j for details), the energy constraint is given by 

Q 



q ' 



(2) 



in which Ti is the hamiltonian of the system under consideration, Si represent the CI possible energy states, and we 
have introduced the concept of escort distribution 



where 



/9 



(3) 



(4) 



being (3 the Lagrange parameter associated with the constraint in Eq.([21), = [1 + (1 — g) x^JJ"^ the exponential 
and [?/]_|_ = y9{y), where 6{y) denotes the Heaviside step function. This implies a cutoff for Si given by 



1-q 



(5) 



The definition of the physical temperature in the nonextensive statistics is still an open issue [25|, |26|, l27|, |28|, |29|, |30|, 
[31I, m, [3^ . From a pragmatic point of view, since has the dimension of energy, (/3^fc)~^ is a temperature scale 

which can be used to interpret experimental results. The validity of this choice was first shown experimentally p^ . 
and later theoretically 0, il3, llB] foi' manganites. 

The MC technique has been successfully used to study the physical properties of Ising models [13, [IH]. Thus, in 
order to generalize the study of the properties of this model by a nonextensive approach, we modified the Metropolis 
method for the nonextensive statistics. To proceed the single spin flip MC calculations [3y| and to obtain the physical 
quantities of interest of the system (magnetization, susceptibility, specific heat, and other quantities), we have changed 
the usual statistical weight to [2ll | 



P, 



i,bef ore 
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(6) 



The above equation is the ratio between the escort probabilities, Eq.Q, before and after the spin flip, and arc the 
energy states related to the Hamiltonian: 



(7) 
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FIG. 1: (Color online) Magnetization versus temperature (left side) and scaling plots (right side) for some values of 0.5 < q < 1.0. 
We observe phase transitions for all values of q in that range at different critical temperatures and with distinct critical exponents, 
that are given in Table I. 



where (ij) denotes the sum over nearest-neighbors on a square lattice of size N ^ L'^, Si = ±1 and J > (which 
imphes a ferromagnetic interaction). 

Since Eq. ([5]) is a ratio, the calculated weight can also be written as the ratio between the g— exponentials with a 
bias q [2l| . It is important to emphasize that Wq is the quantity that will be compared to a random number in the 
Metropolis algorithm and also to note that the cutoff procedure, Eq.(l5]), must be taken into account, i.e., it must be 
included into Wq to avoid complex probabilities. 

Taking into account Eqs. and ([7]), we performed MC simulations with the entropic index q S [0, 1] and linear 

lattice sizes of L = 32, 64, 128, 256 and 512, with periodic boundary conditions and a random initial configuration of 
the spins. The following results were obtained after 10^ MC steps per spin. 
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FIG. 2: (Color online) The pseudo-critical temperatures Tc{L) versus L^^ for some values of q. The extrapolation given us 
the critical temperatures Tc in the thermodynamic limit {L~^ — + 0). Notice that for the q < 0.5 case the Tc{L) values do not 
depend on L. 
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Q 
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0.6 


1.761 ±0.003 


0.34 ±0.01 


0.025 ±0.001 


1.69 ±0.04 


1.00 ±0.01 


0.8 


1.891 ±0.007 


0.15 ±0.02 


0.075 ±0.002 


1.71 ±0.04 


1.00 ±0.01 


1.0 


2.259 ±0.011 


0.00 ±0.00 


0.124 ±0.006 


1.75 ±0.01 


1.00 ±0.02 



TABLE I: Three different entropic indexes in the range 0.5 < q < 1.0 and its respective critical temperatures and exponents 
a, P, 7 and v. For q — 1, the critical exponents and temperature are very close those expected, a — 0.0, /3 — 0.125, v — 1.0, 
7 = 1.75 and Tc — 2.269. Notice that the exponents v are essentially the same for the three cases, whereas a, /3, 7 and Tc 
depends on q. We have found a logarithmic dependence of a on the lattice size L in the q = 1 case, as expected, which give us 
a{q = 1) = 0. The errors in the numerical estimates of the critical temperatures and the critical exponents are also presented. 

III. NUMERICAL RESULTS AND FINITE SIZE SCALING 

In the following we will show results for the magnetization per spin, m, and for the susceptibility x and the specific 
heat C, which can be obtained of the simulations from the fluctuation-dissipation relations, 

(m2)-(m)2 

X = 7^ , (8) 

where ( ) stands for MC averages and e is the energy per spin (we have considered J = fc = 1 for simplicity). In Fig.[T] 
it is shown the simulations in the range 0.5 < q < 1.0, where we present on the left side the magnetization versus the 
temperature for q = 0.6, q = 0.8, and q ~ 1.0. As may be observed on this figure, the magnetization curves changes 
continuously from a ordered ferromagnetic phase to a disordered paramagnetic one. Thus, the critical exponents and 
the critical temperatures of the model can be obtained by the standard finite-size scaling (FSS) forms, 

m{T,L) = L-P'- m{{T-T,)L^'-) , 
X{T,L) = L'^/'' x{{T-T,)L^^n , 
C(T,L) = L"'" C{{T -Tc) L^'") , 

(10) 

where a is a constant. In Eqs. (|10p . the exponent (3 is related to the behavior of the magnetization near the critical 
point Tc, v is related to the divergence of the correlation length, 7 governs the divergence of the susceptibility at the 
critical point, a is related to the divergence of the specific heat at Tj, and m, x and C are scaling functions. The critical 
temperatures of the infinite lattices Tc were obtained by extrapolating the Tc(L) values given by the susceptibility 
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FIG. 3: (Color online) Upper figures: Binder cumulant, Ea. (|lip . versus temperature for q — 0.6 [Fig. 3(a)] and the best collapse 
of data [Fig. 3(b)], based on the FSS in Ea. (|12|l . The parameters are v = 1.0 and Tc = 1.761. Lower figures: magnetization 
values at each pseudo-critical temperature Tc{L) for various lattice sizes L [Fig. 3(c)] and the corresponding values of the 
susceptibility peaks [Fig. 3(d)], for q — 0.6. The fittings in the log-log scale give us the corresponding values of the critical 
exponents ratios P/v and j/v, 0.025 ± 0.001 and 1.69 ± 0.04, respectively. 



peaks positions^ (see Fig. [2]). On the other hand, the exponents f were obtained by means of the Binder cumulant 
[33], defined as 

which has the FSS form 

Ul^Ul{{T-T,)L^''') , (12) 

where Ul is a scaling function that is independent of L. The error bars in the estimations of v were obtained 
following the standard procedure for collapsing data of the Binder cumulant in the finite-size scaling approach, i.e., by 
monitoring small variations around the best collapsing pictures. In Fig[3]we show, as an example, the Binder cumulant 
for q = 0.6 [Fig. 3(a)] and the best collapse of data [Fig. 3(b)], based on Eq.([T2]), obtained with the critical temperature 
Tc and the exponent v given in Table IH Also in Fig[3] we show, for q = 0.6, the values of the magnetization at the 
pseudo-critical points Tc{L) for various lattice sizes L [Fig. 3(c)] and the corresponding values of the susceptibility 
peaks positions [Fig. 3(d)], in the log-log scale. Linear fitting of data yield the parameters: 

13 /v ^ 0.025 ±0.001 , (13) 
j/i^ 1.69 ±0.04 . (14) 

By repeating the fitting procedures of the specific heat peaks positions versus the lattice size L, in the log-log scale, 
we have estimated the ratio a/i' = 0.34 ± 0.01, for q = 0.6 (see Table We have calculated the exponent i' by 



[1] Equivalcntly, we can determine the Tc{L) values by the maxima of the specific heat curves. 
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means of the Binder cumulant, as above-discussed, which allowed us to estimate the critical exponents a, (3 and 7. 
The procedure was the same for the other values of the entropic index and the best collapse of the magnetization 
data, presented on the right side of Fig. [H supports the validity of the FSS forms in Eqs. (fTO|) and the reliability of 
the numerical results for the critical exponents. The obtained numerical results are summarized on Table |T1 Note 
that for (J = 1, the critical exponents a, /?, 7 and v are quite close to the exact known values of the standard 2D Ising 
model, as expected. However, for different values of q in this range, we have found different values of the exponents 
a, (3 and 7, as we can see in Table HI whereas the values of v are the same, within the determined uncertainty. These 
results will be discussed in with more details bellow. 






FIG. 4: (Color online) Results for the magnetic susceptibility x for L = 128 and typical cases of g < 0.5 [g = 0.2, (a)] and 
q > 0.5 [g = 0.8, (b)]. It is also shown the specific heat curves for the same values of g [(c) and (d), respectively]. Although it is 
possible to observe a jump on the left figures, the histograms of the energy states visited during the MC simulation [Figs. 4(e) 
and 4(f)], at the corresponding critical temperatures, show only one-peak structures, indicating continuous phase transitions, 
even for the case q < 0.5. We have defined the energy as the fraction of unhappy bonds in the system, e = {E + 2N) /{AN), 
where E is the total energy given by Eq. ([TJ and A'^ is the total number of spins. 

In Fig. [4] (right side) we show as example the susceptibility and the specific heat for q = 0.8, as well as a histogram 
of the energy states visited during the dynamics of the system, at the critical temperature. This histogram shows 
only one peak, i.e., we have a continuous phase transition, as shown in the magnetization curves, Fig[TJ 

In Fig. Oit is shown the behavior of the magnetization as a function of the temperature for two cases of g < 0.5. 
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FIG. 5: (Color online) Magnetization versus temperature for q = 0.2 and q = 0.5. We observe jumps on the curves at the 
corresponding critical points, but these Tc values are independent of the lattice size, showing that the scaling naturally occurs 
on the q < 0.5 regime. 



One can see that the critical temperatures, Tc ^, are the same for all lattice sizes, as shown in Fig. [2] This result is a 
consequence of the cutoff of the escort distribution, Eq.®, and the magnetization jumps at Tc from m = 1 to m = 
(for more details, see HH). The cutoff also affects the susceptibility and the specific heat, as we can see in Figs. 4(a) 
and 4(c), respectively. However, if we compute histograms of the energy states visited during the dynamics'^, at 
the critical temperatures, we can verify that we have only one peak for all q G [0, 1], indicating the occurrence of 
continuous phase transitions [see Figs. 4(e) and (f)]. In other words, the cutoff keep the MC simulation trapped in 
the ground state for T < Tc and the thermodynamic quantities suddenly change at Tc. 



[2] According to Ref. |21| . the critical temperatures in this regime are given by Tc = 4 (1 — g). 

[3] The energy per spin curves show jumps at the critical temperatures Tc for the q < 0.5 cases, but the histograms of the energy states 
visited clearly show one-peak structures, indicating the occurrence of continuous phase transitions. In other words, the jumps are only 
an effect of the cutoff of the Tsallis distribution, as in the magnetization curves. 
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FIG. 6: (Color online) The critical exponents a, (3 and 7 as functions of the entropic index q. In the range 0.5 < q < 1.0, the 



dependencies on q are given by a{q) = (10 g - 33 q + 23)/20, /3(g) = (2g- l)/8 and 7(g) 
a, P and 7 fittings are less than 5%. 



g + 7)/4. The error bars for 



We can see from Table U that the critical exponent v does not depends on q in the range 0.5 < q < 1.0, and we 
conjecture that the correct value for any g is = 1.0. However, a, /3 and 7 depend on the value of q. Fitting the 
numerical values of a with a second-order polynomial function of we have found that a{q) = 0.5 q" — 1.65 q + 1.15, 
for 0.5 < g < 1.0 (sec Fig.H]), or 

a(9) = ^(10g'-33g + 23), (15) 

which give us the exact known value a{q = 1) = and a[q = 0.8) ~ 0.15 and a{q = 0.6) = 0.34, in agreement with 
the values given in Tab. |T1 In addition, fitting the numerical values of 7 also with a second-order polynomial function 
of g, we have found that 7(g) = 0.25 q^ - 0.25 q + 1.75, for 0.5 < g < 1.0 (see Fig. E]), or 

l{q) = \{q^-q + l), (16) 

which give us the exact known value 7(5 = 1) = |, 7(9 = 0.8) = 1.69 and 7(9 = 0.6) = 1.71. These results are also 
close to the ones obtained numerically. On the other hand, fitting the numerical values of [3 with a straight line, we 
have found that f3{q) = — 0.124 -F0.249q, for 0.5 < g < 1.0 (see Fig. IB]). Wc may conjecture that the exact dependence 
of P on q in this range is 

/3(g) = ^(2g-l), (17) 

which give us the exact known value /3{q = 1) = i, f3{q = 0.8) = 0.075 and f3{q = 0.6) = 0.025, values that are also 
close to ones the obtained numerically. 

These results suggest a nonuniversality of the critical exponents along the ferromagnetic-paramagnetic frontier. In 
addition, it also suggest that the scaling relations 

2/3 + 7 = di^ , 

a + 2/3 + 7 = 2 , (18) 



[4] Wc have found a logarithm dependence of a on the lattice size L in the q = 1 ease, as expected, which give us a(q = 1) = 0. 
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where d is the dimension of the lattice (d = 2 for the square lattice), should be changed. Thus, if we consider the 
above dependence of a, (3 and 7 on g, the first scaling relation of Eqs. p^ . will become 

2/3 + 7 = {d + nq)v , (19) 

where 

= + . (20) 

Notice that for 9 = 1, one has ni = 0, and the standard scaling relation is recovered. On the other hand, although a, 
/3, and 7 depend on the entropic index g, the Rushbrooke equality is satisfied for all 0.5 < q < 1.0, within imccrtainty. 



IV. CONCLUSIONS 



We have studied the Ising model with nearest-neighbors interactions on a square lattice by means of numerical 
Monte Carlo simulations. In our approach, different from other authors [2^, [2^, [H, [s^, , we simply changed the 
weight in the Metropolis algorithm to a ratio between the escort probabilities of the nonextensive statistics. This 
study was motivated by possible connection of the Tsallis statistics and some manganese oxides, called manganites, 
like Lao.eoYo.oyCao.saMnOs [l3, [H, [3 ■ Due to computational cost, our simulations were done after 10^ Monte 
Carlo steps, with the entropic index q G [0, 1] and the linear lattice sizes L = 32, 64, 128, 256 and 512. 

The Monte Carlo simulation of an Ising model with nearest-neighbors interactions showed a distinct behavior of 
the same system considered in the infinite-range-interaction limit [l9| . Jumps on the magnetization and susceptibility 
curves in the range 0.0 < q < 0.5 occur in both approaches, but for short-range interactions we do not have first-order 
phase transitions. In addition, the mean- field calculations foresee the same critical exponents of the 2D Ising model 
in the framework of the Boltzmann-Gibbs statistics. However, our calculation of the magnetization, the susceptibility 
and the specific heat for the short-range interacting system showed that three of the critical exponents depend on q 
in the range 0.5 < q < 1.0. 

Finite-size scaling analysis of the results showed that the critical exponents a, (3 and 7, that are related to the 
behavior of the specific heat, the magnetization and the susceptibility near the critical point Tc, respectively, depend on 
q in the range 0.5 < q < 1.0. Based on the numerical estimates of these exponents, we conclude that the dependencies 
are of the form a{q) = (10 _ 33 g + 23)/20, /5(g) = (2 g - l)/8 and 7(g) (g^ - g -f 7)/4. Although the exponents 
a, (3 and 7 depend on g, as well as the critical temperatures Tc [2l|, the exponent v does not; we found that = 1.0 
V g. These dependencies of the critical exponents on the entropic index suggest a nonuniversality of those exponents 
along the ferromagnetic-paramagnetic frontier. It also suggest a violation of the scaling relations a-|-2/3 + 7 = 2 
(Rushbrooke equality) and 2/3 + 7 — However, when we take into account the g— dependence of the critical 
exponents showed in Tabled we notice that the former scaling relation should be changed to 2 /3 + 7 = (d + Uq) 
where Uq = (g^ + g — 2) /4 (note that for g = 1, we obtain ni =0), but the Rushbrooke equality is not altered. Thus, 
the inhomogeneities introduced in the system by the nonextensive statistics may be responsible for the g— dependence 
of the critical exponents a, /? and 7, as well as the critical temperatures T^. 

On the other hand, we have a completely different scenario in the range 0.0 < g < 0.5. The cutoff of the Tsallis 
distribution keep the system in the ground state (with m = 1) for T < Tc = 4(1 — g), and at Tc the magnetization 
jumps suddenly to zero, i.e., to a equiprobable state [2l] |. In the same way, the susceptibility and the specific heat 
curves also present jumps at Tc, due to the cutoff. Although the presence of these jumps, the histograms of the energy 
states visited during the dynamics, at the critical temperatures, show only one-peak structures, which is a indicative 
of the occurrence of continuous phase transitions. 

Previous works on long- and short-range interactions ID Ising models [4l|, [4^ predict that the magnetization 
scales differently for g < 1.0 and g = 1.0 regimes. Therefore, in this work we showed that the magnetization of the 
short-range 2D Ising model scales also differently in two regimes: for 0.5 < g < 1.0 the system scales as a 2D Ising 
model, but for g < 0.5 the magnetization and the critical temperature are independent of the lattice size due to the 
cutoff; thus, the scaling appears naturally on the system. 

Also in a previous work [25[ , it was shown that 2D Ising model with nearest-neighbors interactions does not undergo 
phase transitions, except for g = 1.0. The main difference between their approach and ours is related to the definition 
of the temperature scale. In that work [25| . the authors have chosen /3 as the parameter related to the temperature 
scale and, in this work, we have chosen /3^. The relation between these parameters is given in Eq.Q. The advantage 
of our approach over previous one, for the choice of the temperature scale, is that ours is supported by previous 
description of the magnetic properties, experimentally and theoretically investigated, of manganites [6l. FtI. [itI. [TsI. [l9j . 
Thus, based on that, we believe that the 2D Ising model undergoes a phase transition even for g 1.0, and the scaling 
relations should be changed as described above. 
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Extensions of this work to describe inhomogencous magnetic systems, i.e., systems in which the exchange interaction 
changes along the sites of the lattice, as well as the study of the effects of uniform and random magnetic fields, within 
a nonextensive approach would be of great interest, because it can yield some clues to questions about the connection 
of such systems and the nonextensive statistics. 
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